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1 Introduction 

The quantum transfer matrix [16, 17] is a means which makes it possible to calculate 
finite temperature properties of one-dimensional quantum spin systems. It becomes 
particularly efficient [7, 8] for so-called integrable models which are related to repre- 
sentations of the Yang-Baxter algebra and are solvable by Bethe ansatz. Among these 
models the spin-^ XXZ chain (henceforth briefly referred to as 'the' XXZ chain) is one 
of the most thoroughly studied (see e.g. [2, 12, 18]). In contrast to many other integrable 
models, not only its spectral properties and thermodynamics are well understood. In 
addition, efficient formulae for norms [11], scalar products [15], and the multiple ac- 
tion of certain non-local operators [5] are available for the generalized model [11] 
defined by the 7?-matrix of the XXZ chain. These technical achievements led to the 
derivation of certain multiple integral representations for (static) two-point correlation 
functions and generating functions of two-point functions at zero temperature [5]. 

In this work we show how the results of [5] can be generalized to finite tempera- 
tures. We concentrate on the generating function of the a^-a^-correlation functions [3]. 
We use the quantum transfer matrix and the formalism of non-linear integral equations 
in order to express it as a multiple integral over the same auxiliary function that deter- 
mines the free energy of the model. In the zero temperature limit we recover the result 
of [5]. An extension of our results to other finite temperature correlation functions of 
the XXZ chain is possible and will be published elsewhere. 

After providing the necessary foundations for treating correlation functions at finite 
temperatures in section 2 we derive our main result in section 3, which is the integral 
formula (116) in theorem 1. As a special case of (116) we obtain an integral repre- 
sentation for the so-called emptiness formation probability [13] of the XXZ chain at 
finite temperature. In order to keep our presentation short we formulate our results in 
the first instance only for the off-critical regime (A > 1) of the model. Note, however, 
that our method applies for arbitrary values of the anisotropy parameter. In section (4) 
which is devoted to an outlook and to conclusions we indicate how our results have to 
be modified in the critical regime. 

2 Foundations 

2.1 Hamiltonian and i?-matrix 

The Hamiltonian of the L-site XXZ chain, 

Hxxz = jt fe-i^ + + A(a}_ia} - 1)) , (1) 

is defined on the L-fold tensor product (C^)®^ through the elementary operators oj, 
a = Xjy,z, which act as Pauli matrices on the jth factor and trivially elsewhere. For 
7 = we set Oq = a", thereby specifying the boundary conditions to be periodic. 
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/ > fixes the energy scale, and the real parameter A is the anisotropy parameter of 
the model. 

The Hamiltonian (1) is closely related to a certain trigonometric solution 
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of the Yang-Baxter equation which not only allows one to diagonalize the Hamiltonian 
by means of the algebraic Bethe ansatz, but will also enable us to associate a quantum 
transfer matrix with the model that will be our principal tool for calculating its finite 
temperature properties. We shall now briefly review the general theory. 

For simplicity we shall restrict ourselves to fundamental models. These are models 
which are completely determined by a numerical d} x matrix RCkjiu) e End(C'^ 
C^) satisfying the Yang-Baxter equation* 

R'^^,MRa^{Ky)RlUM,v) = Rll/f^,v)R'^^,{l,v)R^!^„M . (3) 

This equation is nowadays not usually presented in components, but either graphically 
or in 'operator form' with respect to the gl{d) standard basis^ canonically embedded 
into (C^)^ . In the latter case one defines 



RjkM=RllMejlek'y. 

Then (3) is equivalent to 

RnMRl3MR23{lil,v) = R23{lil,v)Ri3MRi2{X,lil) . 



(4) 



(5) 



A Hamiltonian associated with the 7?-matrix R{X,jii) is usually constructed as fol- 
lows: First of all an L-matrix at site j with matrix elements 



Lj^M = R;iMej', 

is introduced. These matrix elements are operators in (End(C^))®^. 
the Yang-Baxter equation (3) by ej'!^' implies that 

RMiLjil,w)^Lj {^,v)) = {Lji^,v)0Ljil,v))RM, 

*Here and in the following we sum over doubly occurring Greek indices. 

^The basis consists of all <i x <i matrices having a single 1 at the intersection of the ath row with 
the pth column. Accordingly the basis matrices multiply as e^ey = S^e^. The canonical embedding is 
defined by ej^ 



(6) 
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j) ^-j-j^ being the <i x J unit matrix. 
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where = i?pg. One says thatLj(A,, v) is a representation of the Yang-Baxter algebra 
with i?-matrix i?(A,,//). This representation is called the fundamental representation. 
Next, a monodromy matrix T{X) is introduced, 

r(X)=Li(X,v)...Li(X,v). (8) 

Since [Lj+i^{X,v),Ljl{X,v)] = it follows that 

R{l,^){T{l)®T{^i)) = {T{p)^Til))R{l,^). (9) 

This means that also T{X) is a representation of the Yang-Baxter algebra, a property 
which we use synonymously with integrability. 

We shall assume that is regular in the sense that ^(0, 0) = P is the permuta- 

tion matrix acting on (8) C^. Then it follows that Lj'^{0, 0) = e^p and further, setting 
V = in (8), that the transfer matrix, defined as 

?(X) = tr(r(X)), (10) 

has the expansion 

t{'k) = Uexp{'kH + 0{}?)} (11) 
around A, = 0. Here U = t{0) is the shift operator and 

(12) 



7=1 



A,=0 



is the Hamiltonian associated with the fundamental representation (6) of the Yang- 
Baxter algebra. Note that periodic boundary conditions, ^oi = ^Li , are implied in ( 12). 

Employing the usual conventions for tensor products we may interpret the matrix 
(2) as a linear operator in C^. It is then easy to show that it satisfies the Yang- 
Baxter equation (3). Inserting (2) into the right hand side of (12) we obtain 



^ = 2^ ^ (a5_ia^ + a^._ia^ + ch(Ti)(a^._,a^- 1)) . 



(13) 



sh(ri) j 
Thus, 

Hxxz^2Jsh{n)H (14) 

if we identify A = ch(r|). 

2.2 The Trotter-Suzuki formula 

Consider the elementary formula 

lim(l + f)''=e^, (15) 
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valid for any sequence (Xiv)iveN of complex numbers converging to X. It still works 
[16] for a matrix (an operator) X = —^H. Hence, (11) and (15) imply that 

taJtZ-'K-l))" = + 1 (-p// + (1)))" = e-P" . (16) 

This is a special form of the Trotter- Suzuki formula for the statistical operator 
if we choose [3 to be the inverse temperature, P = l/T. 



2.3 The quantum transfer matrix 

We would like to use the Trotter- Suzuki formula (16) in a way that is compatible with 
the special structure of integrable models. In this context it will turn out to be useful to 
introduce A'^ auxiliary spaces denoted by indices l,...,N and to work with i?-operators 
rather than L-matrices. 

We introduce two types of monodromy matrices: type one, 

Tj{1)=RjlM...RjiM, (17) 

and type two. 



Tjil)=Rij{,u,l)...RLjM. (18) 



Remark. 



If ^ is unitary, ^i2(?i,^)i?2i(j",?t) = id, then T j{X) = T-7'{X) 



By construction the trace of the monodromy matrix Tj{X) with respect to the aux- 
iliary space is the usual transfer matrix, t{X) = tr j{Tj(k)). The structure of f(X) = 
tvj{Tj{'k)) can most easily be explored by means of the parity operator P which we 
define by its action on the canonical basis of local operators, Pef^P = e^-y+ip. Then 

Tj{l) = P{PRP) j^i^.l) . . . {PRP)j,{^,l)P. (19) 

Here one has to be careful not to confuse the parity operator P with the permutation 
matrix P. In (19) we set = 0, take the trace in space j and take into account that the 
parity operator acts only in quantum space. We infer with (1 1), (12) that 

= PUcxp{l £ d^{PRP)j_ij{0,l)\^^^ + ...}P 

= exp{X £ d^Rj.ij{0,l) l^^o + 0{1^)}U-' , (20) 

where we used in the second equation that {PRP)j-ij — Rjj-i and PUP = U~^. 
In the following we shall assume that 



(21) 
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This is certainly true for /?-matrices of difference form like (2) but also for other more 
general i?-matrices like the i?-matrix of the Hubbard model. Combining (11), (12), 
(20) and (21) we obtain 



P».L := (<(l)'(-|)) ' = (l + K-P'' + " (^))) ' ■ <22> 
Using the second equation (16) we conclude that 

limp^,z. = e-P^. (23) 

iV— >oo 

The latter expression is very useful for integrable models. It leads us rather directly 
to the notion of the quantum transfer matrix, for we may write 

PA^,L = tri...^{T^(|) %3t(-|) . . . T, (-1) } 

= tri...47^Hl)^fcT(-|)-^2(|)7'f(-|)} 
= tri...yv{^iyv(A', |)^2yv(A', |) • • •^Lyv(A', |) 

^12(i"J)^22(i"J)---^L2(i",|) 



A/=0 



= tv-,_,_^{T,^™{0)...T^^''{0)}. (24) 

Here the transpose with respect to space 1 occurring in the third equation is defined 
by R^^^l = R^- Notice that we have reordered the product of i?-matrices in the fourth 
equation and have introduced another monodromy matrix, 

7;P™W=^,s(X,|)fl;L^,(-lX)...fl,,(X.l)«i'(-iX) (25) 
in the fifth equation. This monodromy matrix and the corresponding transfer matrix 

tQ™{X)=trjT^^''{l), (26) 

which will be called the quantum transfer matrix, are our main tools for exploring the 
finite temperature properties of the XXZ chain below. 
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The quantum transfer matrix is useful mainly for two reasons. First, in the ther- 
modynamic limit L — > oo a single leading eigenvalue of the quantum transfer matrix 
determines the free energy. Second, the spectrum of the quantum transfer matrix at 
finite Trotter number A'' can be calculated by means of the algebraic Bethe ansatz. Let 
us explain this in more detail: 

Using (24) the partition function for a quantum chain of length L is calculated as 

oo 

Zl = lim tri..xPyv,L = Hm tri...;v(^^™(0))'' = £ A^(0) , (27) 

where the A„(0) are the eigenvalues of the quantum transfer matrix in the Trotter limit 
— > oo at spectral parameter A, = 0. The free energy per lattice site in the thermody- 
namic limit follows as 

if Ao(0) is the (finitely degenerate) leading eigenvalue of the quantum transfer matrix. 
For the remainder we shall accept the following. 

(i) The limits L ^ oo and — > oo are interchangeable [16, 17]. 

(ii) The leading eigenvalue of the quantum transfer matrix is non-degenerate, real 
and positive and in the Trotter limit is separated from the next-to-leading eigen- 
values by a gap. 

The thermodynamics of the model under scrutiny is then determined by the leading 
eigenvalue Ao(0) of the quantum transfer matrix in the Trotter limit N ^ oo. 

We have seen that instead of calculating the infinitely many eigenvalues of the 
usual transfer matrix which determine the spectrum of the Hamiltonian it suffices to 
calculate a single eigenvalue of the quantum transfer matrix. Now we shall see that 
the quantum transfer matrix is a representation of the same Yang-Baxter algebra as the 
ordinary transfer matrix. Taking the transpose with respect to space 1 in (5) we obtain 

R23{:K,i^)R'l2{v^)R'l^{v,i^)=R'l^{v,i^)R'l2{v,X^^^^ (29) 

Then, because of (5) and (29), the monodromy matrix (25) provides a representation 
of the Yang-Baxter algebra, 

i?,,(A,^)7;e-^(A)rf ^^(^) = rf^^(^)r,e^^(A)^ .,(A,^) . (30) 

This means that the leading eigenvalue of the quantum transfer matrix may be obtained 
by algebraic Bethe ansatz (or possibly by one of the other powerful methods connected 
with the Yang-Baxter algebra as e.g. the method of separation of variables [14]). 
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2.4 Correlation functions within the quantum transfer matrix ap- 
proach 

The quantum transfer matrix approach allows us to calculate, in principle, finite tem- 
perature correlation functions of local operators by means of the algebraic Bethe ansatz. 
Here we shall consider correlation functions of the form 

{xl'K..xt'^\=nn. "' '\ , (31) 

tri,...,Le~P^ 

(n) 

where j,k e {I, . . .,L}, j <k and the are arbitrary local operators. Using (23) and 
(24) we rewrite (31) as 

(xf\..xr^\ 

= lim tr,,^^str,.x7-«™(0)...Te™(0)xj»...xf-'+'Vz„.. 

Jy — ^oo -f 

= lim trj ^(rer^(0))^"'tr{rer«(0)x(i)}...tr{re^^(0)x(^->+i)} 

N,L—K=° 

L^!o^An~^+'"~\o)(yn|tr{r6™(o)x(i)}...tr{rQ™(o)x(^-^+i)}|y„) 

= lim A^"^"^(0)(^o|tr{r2^^(0)z(i)}...tr{r2™(0)x(*-^'+i)}|^o)- (32) 

Here we assumed that the quantum transfer matrix is similar to a diagonal matrix* and 
has 'normalized' eigenvectors |^n). We see that a single normalized eigenvector |^o) 
and the corresponding eigenvalue Ao{X) of the quantum transfer matrix completely 
determine the finite temperature correlation functions of local operators. 

Let us consider two examples related to the XXZ chain. In this case d — 2, and we 
can represent the monodromy matrix (25) as a 2 x 2 matrix 

^''™(^) = (c(X) DW) '''' 

in 'auxiliary space'. 

Example. X^^) = a , = o+, and all other X^") are equal to the 2 x 2 unit 

matrix I2. Then (32) turns into 

(070+)^ = lim Ar'"'(0) (^o|5(0)(A(0) +D(0))^-^-iC(0)|^o) • (34) 



*This is known for the XXZ chain, but not in the general case. 
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Example. = z(2) = • • • = = e^^2 = ( J o, ) and 7 = 1, A: = m. In this case 

(32) turns into 

<exp{(pL-=i^'„2|^^ ^ ^iiJL^'^W (^o|(A(0) + e«PD(0)n^o) • (35) 

This special correlation function [3] is a generating function of the c^-o^ correlation 
functions which can be expressed as 



(36) 

(p=0 



where is the 'lattice derivative' defined on any complex sequence (fl„),js=N by 
DmCim — am — cim-\- The function (35) generates a number of other interesting cor- 
relation functions that will be discussed below after we have derived an integral repre- 
sentation for (35). 

The main goal of this work is to represent the generating function (35) as a multiple 
integral. The reason why we begin our study of finite temperature correlation functions 

with the generating function (35) and not with any two-point function like (34) is that 
(35) is slightly more simple than the other interesting correlation functions, because 
the individual factors in the product of operators on the right hand side of (35) are 
closely related to the quantum transfer matrix. 

The expectation value on the right hand side of (35) can not easily be calculated 
directly. Instead, we shall introduce an inhomogeneous version of it and shall perform 
the homogeneous limit only at the end of our calculation. We introduce a 'twisted 
quantum transfer matrix' of the XXZ chain as 

t^{X)=A{X) + Q'^D{X). (37) 

For brevity we have suppressed the label QTM here. It follows from the Yang-Baxter 
algebra (30) that the twisted quantum transfer matrices ?(p(A,) for fixed (p form a com- 
mutative family, 

[t^{X),t^{^l)]=0. (38) 
For any set {^} = {^j}J=i, € C, we define a function 

m m 

^Ni^m) = (^oi [u^m [U'o\^j)] 1^0) (39) 

depending symmetrically (because of (38)) on the ^y. Then the generating function 
(35) can be calculated as 

<exp{(pL-=ie„i}>^= lim lim 4>^(cp|{^}). (40) 

Our strategy below will be to first calculate ^>;v(9|{^}) by means of the algebraic 
Bethe ansatz and then transform the resulting expression into a form that will allow us 
to perform the homogeneous limit and the Trotter limit on the right hand side of (40). 
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2.5 Algebraic Bethe ansatz for the quantum transfer matrix 

We shall now recall how to calculate spectrum and eigenvectors of the quantum trans- 
fer matrix by means of the algebraic Bethe ansatz. The result of an algebraic Bethe 
ansatz calculation is usually not explicit but expresses eigenvalues and eigenvectors in 
terms of the solutions of a set of Bethe ansatz equations. To actually solve the Bethe 
ansatz equations is another problem we deal with in the next section. 

For the algebraic Bethe ansatz the form (9) of the Yang-Baxter algebra is more 
convenient than (30). Hence, we shall write the monodromy matrix (25) of the quan- 
tum transfer matrix as a 2 x 2 matrix. For this purpose we first of all introduce two 
L-matrices L and L setting 

eaUj^M = RajM , eaUjli-n,^) = R%{-fi,'k) . (41) 

Comparing these equations with the explicit expression (2) for the i?-matrix of the 
XXZ chain we obtain the explicit forms of the two L-matrices, 

where we have employed the shorthand notations 

^ '^^ sh(X-// + T|)' ^ '^^ sh(X-// + T|) 

When expressed in terms of the L-matrices (42) the monodromy matrix of the quantum 
transfer matrix becomes 

tQ™{X) = Ln{X, |)L;v-i . . .Li . (44) 

This is now a 2 x 2 matrix of the form (33) in 'auxiliary space' which satisfies the 
defining relations of the Yang-Baxter algebra of the form (9). Moreover, T^^'^iX) 
acts as an upper triangular matrix on the vector 

|0) = (ei®e2)«T = (J)(^(0)(^(J)^...(0) (45) 

^ V ' 

factors 

which is obvious from the form of the L-matrices (42). More precisely, 

C(A.)|0) = 0, A(X)|0) = fl(X)|0), D{X)\Qi) = d{X)\Q) . (46) 



where 



"W=M-|A)^=(^^^f, (47a, 
Vsh(X+^-Tl)/ 

.(X) = Ma)*-(-^^^f. (47b) 
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The first equation (46) is sufficient for the algebraic Bethe ansatz to work. It can be 
performed for arbitrary pseudo vacuum eigenvalues a(k) and d{X} of the diagonal ele- 
ments of the monodromy matrix. These eigenvalues are often called 'the parameters of 
the algebraic Bethe ansatz'. Here they are given by (47). The parameters of the Bethe 
ansatz completely specify the solution which for the general case can, for instance, be 
read up in [12]: 
The vector 

im) = l{Mf=i)=5(^i)---5(M|0) (48) 

is an eigenvector of the quantum transfer matrix ?o(^) if the rapidities (or Bethe roots) 
Xj, j = 1, . . . ,M, satisfy the system 

a(?l,)_-^ sh(X,-X, + Ti) 

d{Xj) |Jsh(X,--X,-Ti) ^ ^ 

of Bethe ansatz equations. The corresponding eigenvalue is 

^^^P^^UX-Xj-^) ,,,-^sh(X-X, + ri) 

Note that the 'creation operators' B{Xj) mutually commute. Therefore the eigenvector 
(48) is symmetric in the rapidities, whence our notation | {X}). 

2.6 Leading eigenvalue and auxiliary function 

In order to calculate the approximant 4>7v(cp|{y) to the generating function (35) of 
the a^-a^ correlation functions we need to know the specific solution of the 

Bethe ansatz equations (49) that determines the leading eigenvalue Ao(X,) and the cor- 
responding eigenvector |^o)- This solution cannot be obtained in a closed analytic 
form. Even worse, unlike in the case of the ordinary transfer matrix, the distribution of 
the Bethe roots cannot be approximated by a smooth 'density function' in the Trotter 
limit, since they accumulate in the vicinity of zero as increases. A way out of these 
difficulties was proposed in [7, 8]: One associates an auxiliary function 

J(X)^ ^h(X-X,-FTl) 
a(A)^^-\sh(X-X^-Ti) 

with the solution {A,} = {X^}^^^ of the Bethe ansatz equations that determines the 
leading eigenvalue of the quantum transfer matrix. It then turns out that the function 
(x{X) is sufficiently well determined by the gross properties of {X}. Knowing that the 
leading eigenvalue for fixed Trotter number A'^ is determined by 

(i) M = N/2 Bethe roots§ which 



^This was already used in the definition (51) of a(^). For technical reasons we assume that s 1 
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(ii) are mutually distinct and 

(iii) lie all on the imaginary axis 

we can derive a non-linear integral equation which (alternatively) defines a(?i). The 
details of the derivation have been described on several earlier occasions [9]. Here 
we only note that (i)-(iii) can be established 'perturbatively' in the high temperature 
limit and have been verified numerically for a wide range of Trotter numbers and 
temperatures. Conditions (i)-(iii) determine the analytical properties of the auxiliary 
function o(X,). 

(i) The auxiliary function a is meromorphic and periodic with period in in imagi- 
nary direction. 

(ii) It has 3N/2 poles (including multiplicities) in the fundamental domain D = 

{z e C| ~k/2 < \lmz\ < n/2}: N/2 simple poles at Ij+t], j = I,..., N/2; 
an A^/2-fold pole at = -^/N; and an A^/2-fold pole atX = ^/N-^. 

(iii) The meromorphic function 1 + a has 3N/2 zeros in D: N/2 simple zeros at the 
Bethe ansatz roots, 

l + a{lj) = 0, (52) 

clustering close to X = as |3 ^ 0; A'/2 zeros which are close to T] as P ^ 0; 
and N/2 zeros which accumulate at — r| for (3 — 0. We assume of the latter two 
kinds of zeros that the modulus of their real part stays larger than Ti/2 for all 
finite temperatures. 

Once the analytical properties of a and 1 + a are established it is not difficult to 
derive the following integral equation for the auxiliary function. 



In a{X) — In 



sh(X-|)sh(^+^+Tl)''' 



,sh(X+|)sh(X-|+Tl), 



^ f d& sh(2ri)ln(l + a((D)) 
c 2ni sh(A,-a)-|-ri)sh(A,-co-ri) ' 



where (for A > 1) the contour C (see figure 1) is a rectangular contour with edges 
parallel to the real axis at ±i7r/2 and to the imaginary axis at ±Y where < y < ri/2. 
Equation (53) is valid inside the rectangle Q defined by |Re?i| < r\/2 and |Im?i| < n/2. 
This means that it defines the auxiliary function a(k) for all other Xin Q once it has 
been calculated on C. In particular, one may calculate the zeros of 1 + o inside C 
which, by construction (see (52)), are those solutions of the Bethe ansatz equations 
that determine the leading eigenvalue and the corresponding eigenvector. We will, 
however, not use equation (53) in this way. It turns out that all quantities we are inter- 
ested in, the free energy and the generating function of the a^-a^ correlation functions 
in particular, can be expressed as integrals over the auxiliary function a on the con- 
tour C. 

It is very important that the contour C of the integral on the right hand side of (53) 
does not depend on N. The Trotter number N merely appears as a parameter in the 
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Figure 1: The canonical contour C for the off -critical regime A > 1. 



inhomogeneity of the integral equation. Hence, the Trotter limit can be performed in 
(53). We obtain 

l^niX) = 27sh2(r|) _ I d& sh(2Ti)ln(l + a(a))) 

^ ' rsh(X)sh(x+'n) ic 27ti sh(X-to+'n)sh(X-to-'n) ■ ^ ^ 

Here we already took into account the rescaling of the Hamiltonian (see (14) and (16)) 
and introduced the temperature T as 

r = H^. (55) 

The auxiliary function a encodes the information about the location of those Bethe 
roots that determine the leading eigenvalue of the quantum transfer matrix and the 
corresponding eigenvector. Using the results of sections 2.3 and 2.4 we conclude that 
a thus encodes not only the complete information about the free energy but also all 
the 'thermodynamically relevant information' about the correlation functions of local 
operators. 

2.7 Including the magnetic field 

Before coming to the calculation of physical quantities in the next sections we have to 
discuss the issue of how to include a (longitudinal) magnetic field into the formalism. 
For this purposes we have to replace the statistical operator e~P^ by , 
where h denotes the magnetic field and 

S' = \Y^<5] (56) 
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is the operator of the z-component of the total spin. 

The reason why this is a rather simple task lies in the commutativity of the i?-matrix 
(2) with the tensorial square of the diagonal matrix = diag(e^,e~®), 

[R{X,lu),@^@] = 0. (57) 

Taking the derivative of this equation with respect to at 9 = we obtain 

[i?(;i,//),00/2+/2 0] =0. (58) 

The latter equation implies that the ordinary transfer matrix and hence the Hamiltonian 
commutes with S^. On the other hand (57) is equivalent to 

which means that is a spectral parameter independent representation of the Yang- 
Baxter algebra with i?-matrix (2). It follows that 

limpArLe^=e-(P«-^^Vr)^ (60) 

while, on the other hand, 

PN.Le^ = 7-e™(0) , . . .2f ™(0) («7 ^_»,, (61) 

where rp^^(X) is the monodromy matrix corresponding to the quantum transfer ma- 
trix as defined for zero magnetic field. 

Thus, we can include a longitudinal magnetic field by simply replacing 

(62) 

in our former calculations. The important point here is that, because of (59), the re- 
defined monodromy matrix is still a representation of the Yang-Baxter algebra with 
the same i?-matrix (2). Therefore our former calculations go through without almost 
any modification. We simply replace a(X) by a(?i)e''/^^ and by (i(?i)e^''/^^ in 
(49)-(51). The only point one then has to take care of is that the loci of the Bethe roots 
that describe the leading eigenvalue are now off the imaginary axis. As a consequence 
the parameter y appearing in our definition of the path above has to be chosen large 
enough (but still between zero and T|/2). Then the only modification of the integral 
equations (53) and (54) is the appearance of an extra term —hjT on the right hand 
side. For example, (54) turns into 

l„^m^ ^ 27sh2(ri) /■ sh(2ri)ln(l+a(co)) 

^ ' T rsh(;i)sh(;i+'n) yc27iish(X-to+'n)sh(X-to-Ti)' ^ ' 

This is the non-linear integral equation for the auxiliary function o of the XXZ chain 
exposed to a finite, longitudinal magnetic field h. 
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It is often useful (in particular when it comes to numerical calculations) to con- 
sider the auxiliary function o together with its 'dual' a = 1/a. As we can see from 
equation (52) the function a is equally useful in determining the Bethe roots, for 
1 + a(k) = is equivalent to 1 + a(?i) = 0. An integral equation for a can be di- 
rectly obtained from (63). For this purpose we choose X and to from inside the 
rectangle Q. Then |ReA,|, |Reto| < r\/2 and hence A, — to±r| 7^ 0. Taking this into 
account we change the free variable in equation (63) from X to CO, multiply it by 
sh(2r|)/27ii sh(A, — CO + r|) sh(A, — co — r|) and integrate over C. The resulting formula 
is 

Jco sh(2Ti)ln(a(co)) 2/sh(T|)sh(2Ti) 



'c 27ii sh(A,-co + T|)sh(A,-co-T|) rsh(A,+T|)sh(A,-T|) " 
We use it in (63) and obtain an integral equation for a: 



(64) 



Innri) = ^ 27sh2(ri) r d(o sh(2ri)ln(l + a(cD)) 

T Tsh{X)sh{l-r\) Jc2nish{X-ai + r\)sh{X-ai-r\)' 

Comparing this with (63) we see that the formal difference consists in replacing h and 
r| by —h and — r|. 



2.8 Free energy and magnetization 

Starting from the expression (50) (properly modified such as to include the magnetic 
field) we can derive expressions that represent the free energy per lattice site / as 
integrals either over or a (for details see [9]), 



f(h T) = ---Tf sh(ri)ln(l + a(co)) f ^to sh(ri)ln(l + a(co)) 

' ' 2 JrlTii sh(a))sh(co + ri) 2 Jrlni 



c2%i sh(a)) sh(co + ri) 2 Jc2%i sh(co) sh(co-ri) 



• (66) 



From these equations the thermodynamic properties of the XXZ chain can be cal- 
culated. For example, the magnetization per unit length of the chain at finite magnetic 
field h and finite temperature T, defined as 

m{h, T) ^ \im (67) 

can be calculated by means of the thermodynamic relation 

m{h, T) = ^ — . (68) 



Introducing the function 
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and taking the derivative of (66) with respect to the magnetic field we obtain 

J(o sh(r|) o((o) 



m{hJ) = -\-j 



Ini sh(co) sh(co - rj) 1 + a(co) 



(70) 



Instead of calculating the explicit /z-dependence of the functions o or o and from this 
the function o, we may calculate it from the linear integral equation 

am = ]+f— ^h(2ri) o{(0) 

27li sh(X - to + Tl)sh(X- to -11)1 + 0(0)) ^ ^ 

obtained from (65) by taking the derivative with respect to h. Equation (71) is numer- 
ically better conditioned than (69) and therefore particularly useful if one is interested 
in numbers. 

For later use we derive another formula for the magnetization using an analogy 
with the dressed function formalism for the density functions at zero temperature. We 
define a function G{X) as the solution of the linear integral equation 

G(X) = ^h(ri) f d(D sh(2ri) G{(d) 

sh(X)sh(X-ri) 27ii sh(X-(0 + ri)sh(X-(0-ri) l + a((o) ■ ^ ^ 



Then it is not difficult to show that 

<i(0 sh(ri) o(ca) f d& G((o) 



Ic 2%i sh(to) sh(to - T|) 1 + o(to) Jc 27ii 1 + o(to) 
and it follows from (70), (72) that 



(73) 



1 /■ d(£) G((o) 1 f d(£) G((o) 

m(h, T) = / = -+ • (74) 

^' ^ 2 7^2711 l + a(to) 2 7^2711 l + a(to) 

3 Results 

3.1 Combinatorics 

We shall now derive our main result which is an integral formula for the generating 
function (35) of the o^-a^ correlation functions of the XXZ chain. The general strategy 
we are going to pursue is to first calculate the inhomogeneous, finite Trotter number 
approximant 4>A?((p| {^}) , equation (39), and then transform it into a form such that the 
limits in (40) can be performed. 

The Bethe ansatz eigenvectors (48) are not 'normalized'. Hence, 4>jv((p|{^}), when 
expressed through the Bethe ansatz eigenvectors takes the form 

f i^t'v.) (Wl[n7^i'^fe)][nT.i'o'(4;)]|{^) 

*'''*|'^»- mrn) ' 
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where the in {k} are those that characterize the leading eigenvalue. The dual 

to \{k}) can be expressed through the action of the operators C{X) (see (33)) on the 
dual pseudo vacuum (0| which is the transposed of |0), 

({X}| = (0|C(Xi)...C(X^/2). (76) 

In order to evaluate the numerator in (75) we calculate in the first instance the left 
action of 117=1 ^<p(^j) this task we can largely follow [5]. The basic idea 

is to use the commutation relations 

C{^)A{X) = f{fi,'k)A{X)C{iu) -g{iJ,X)A{n)C{X) , (77a) 
C{ij)D{K) = fMD{K)C{ij) - gMD{n)C{X) (77b) 

contained in the defining relations (9) of the Yang-Baxter algebra in order to move the 
operators A and D (implicit in the expression t(p) in 

(0|C(;il) . . .Cil^/lM^l) ■ ■ M^rn) (78) 

to the left and replace them with their pseudo vacuum expectation values a and d 
(which are the same for right action on the pseudo vacuum and left action on its dual 
(see (42))). Our notation employed in (77) is 

sh(X-;/ + ri) sh(ri) 
^(^'^)= sh(X-^) ' ^(^'^) = ^h(X^- ^'^^ 

Except for the commutation relations (77) one needs in the actual calculation the mu- 
tual commutativity of the families of operators C{X) and t(p{X), respectively, 

[C{X),C{ij)] = [t^{X),t^{lj)] = 0, (80) 

which also follows from (9). 

When we apply (77) successively in (78) in order to move the operators A and D to 
the very left, then the number of operators C per term in a single application of (77) is 
conserved. Merely the arguments of the commuted operators are exchanged in every 
possible way. To state this more formally we need to introduce several set theoretical 
notations. 

We shall deal with finite sets of complex numbers such as For such sets 

we already introduced the short hand notation {X} = {Xyjjlj. Similarly we shall 
write{^} = etc. The number of elements in these sets will be somewhat 

loosely denoted as |^|, ...We further need partitions of sets into disjoint sub- 
sets. By pzi"^} we shall denote the set of all ordered pairs of subsets 
{1+}, {I-} C {1} which satisfy U {^-} = {1} and n {^i^} = 0. The ele- 
ments of the subsets {A,^} will be denoted Xj", such that we can write {A,^} = {A,y^}|^iLi 
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if we want to emphasize the number of elements in the subsets. Using our new notation 
we conclude with (77) and (80) that 



({Min^f(^y)= I Rim\{^-m+}\{i-}){mu{i-}\ m 

j=l {{l+},{X-})ep2{X} 
\^+\+\X-\=N/2 

with numerical coefficients R depending on the sets {^^}, {^^}- These coefficients 
can be calculated in the general situation of arbitrary vacuum expectation values a, d 
of the diagonal elements of the monodromy matrix and for arbitrary sets of mutually 
distinct complex numbers {^}, {A,}, where the do not necessarily satisfy the Bethe 
ansatz equations. 

Lemma 1. [5] The coefficients in (81) are given by the 

formula 



k=l l=\ 



, (82) 



where the so-called highest coefficient 5'({^+}|{A,+}|{A, }) can be expressed as the 
ratio of two determinants: 



s{m\{k-^m-}) = 



detM(X+ ^+) 
detV(X+^+) 



(83) 



with 



^^^^•'^^^-sh(^+-Xt 



l=\ m=l 



7=1 



m=l 



sh(ri) 



sh(;i-^)sh(;i-^+Ti)' 



(84) 



(85) 
(86) 



19 



We may now insert (81) into (75) and use the fact, that the vector |{X,}) in (75) is 
the eigenvector of the quantum transfer matrix that belongs to the leading eigenvalue 

N/2 N/2 

Ao {X) = a{X)Y[ f{Xj, X)+d {X) n /(X, Xj) (87) 
;=i 7=1 

(see (50)). We obtain 

R(ih+\\ih-\\n+\\n-\Mih+\[\i'kr\\i'k\\ 

(88) 



Rmm-}\{x-^}\{x-}){mi}{x-}\{x}) 




niaiAo(^;)^ 





^iv(9|{^})= I 

{{k+}A^-])ep2{X} 

\%+\+\'K-\=N/2 

and only the ratio ({^^} U {^^}|{^})/({^}|{^}) remains to be calculated. But here 
another combinatorial result for the XXZ chain applies. 

Lemma 2. [15] Let {Xyjj^^j a solution of the Bethe ansatz equations (49) and {/Jk}f=i 
any set of distinct complex numbers. Then 



dttN{Xj,^ik), (89) 



{0\C{fii)...C{fiM)B{Xi)...B{XM)\0) 

Ui<j<k<MM'^j - W Mi^k-i^j) 

where 

N{Xj,^,) = K^;.^.) -K^.,^;0STn77f^ • (90) 

^{m i=ij{Xi,Mk) 

In order to apply lemma 2 to our task of calculating ({^"'"} U {X~}\ {X}) / {{X} \ {X}) 
we first insert the set of Bethe roots that characterize the leading eigenvalue of the 
quantum transfer matrix into (89), (90). Then M = N/2, and the right hand side of 
(90) can be expressed in terms of the auxiliary function (51), 

N{Xj,iUk) = t{Xj,Hk)-t{nk,Xj)a{nk) ■ (91) 

Taking account of the Bethe ansatz equations a{Xk) = — 1 it is easy to study the limit 
A'/t — ^ Xk in (89): The prefactor is regular in this limit. Concerning the matrix TV only 
the k\h column is affected, 

lim NiX,,k) = K"^+ , , . . (92) 

n-.X, ^ ' " a{Xk) ^hiXj-Xk + yxjshiXj-Xk-'ix) 

Adopting for a while the notation 

for7=l,...,|^+| 

for 1^+1 + l,...,A^/2, (93) 
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we infer from (89)-(92) that 

{{^+}u{x-}m) 



mm) 



1=1 a{Xj) 



.j=n+lk=l fi'^jj'^k) 



shCkj — Xk) 



<j<k<n sh(^y - ^k) 



detNn 
detA^O 



(94) 



where 



+ ■ 



sh(2ri) 



k= I,. ..,n 
k = n+l,...,N/2. 



(95) 



aCkk) sh(Xy -Xk + r[) sh(Xy - X^t - r]) 

The observation that the columns of N„ and A^o are identical for k > n allows one 
to simplify the ratio of the two determinants on the right hand side of (94). In the 
appendix we prove the following 

Lemma 3. The ratio of the two determinants on the right hand side of (94) is propor- 
tional to the determinant ofannxn matrix, 



detNn 
detNo 



detG(X+^+), 



where the function G(A,,^) is the solution of the linear integral equation 

dG) sh(2ri) G{()),Q 



Gil,Q^ti^,l) + 1 



(96) 



(97) 



Ic 2%i sh(X-a) + ri)sh(X-a)-Ti) l + o(co) 

to be solved on the same canonical contour C as the non-linear integral equation for 

the auxiliary function a. 

Notice that G(?i,^) is an inhomogeneous generalization of the 'finite temperature 
density function' G{X), equation (72), that determines the magnetization (74). 

Using the lemma in (94) and switching back to our original notation we obtain the 
following result. 



({^+}u{X-}|{M) 



mm) 



^ a(ffl(l + a(ffl) 



1^ sh(^;-^/:+Ti) 



n 



.y=U=l/(^7'^fe) 

sh(x+-^+: 



detG(Xt,^+). (98) 



Finally, we use lemma 1 and equations (87), (98) in (88) and remove the explicit 
dependence on the Xj from the resulting expression. The latter is possible by proper 
use of the auxiliary function a and the Bethe ansatz equations (49). The calculation is 
elementary but slightly tedious. We obtain 
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Lemma 4. The finite Trotter number approximant ^>Af((p|{^}), equation (75), has the 
following representation as a sum over partitions of the Bethe roots {X}, that charac- 
terize the leading eigenvalue, and of the inhomogeneity parameters {^}, 



({X+},{X-})ep2W 

|^+| + |X-|=iV/2 



%i({^+}i{^n)2i^+i({^+}i{^+}i{^-}) 



nHa'(xt) nS (1+^(^7)) 



(99) 



where the two functions F„({A,''"}|{^''"}) flnJZ„({A,'^}|{^+}|{^ }) are defined as fol- 
lows: 



n 



sh(Xt-^++Tl)sh(^+-^+-Tl) 



with 



m 1 



,iiish(^;-^++Ti)sh(^;-x+-Ti). 

detM(;it,^+) detG(;i+ ^+) (100) 

(101) 



~ ^ ^ ^ ^ ^ ^ Ashat-X+-ri)sh(X+-^++ri) 



anJ G(A,,^) J5 solution of the linear integral equation (97). 



f{k]Xt)f{K^^] 



), 



(103) 



3.2 From sums to integrals 

In the final step of our calculation we would now like to express the sums over parti- 
tions in (99) as multiple integral over certain canonical paths. We shall apply a similar 
rationale as in [5,6]. 

Consider a function /(tOi , . . . , tOn) : C" — > C, symmetric in its arguments, equal to 
zero if any two of its arguments agree, and analytic on and inside the simple n-fold 
contour C". The function 1 + o(€o) is meromorphic inside the rectangle Q (see section 
2.6), where its only zeros are the simple zeros located at those Bethe roots which 
characterize the leading eigenvalue of the quantum transfer matrix. Hence, the only 
poles of the meromorphic function 1/(1 + ci(to)) inside C are simple poles at the Bethe 
roots with residues 



1 



(to)i la)=A,^- a'(A,y) ' 



; = l,...,A^/2. 



(104) 



22 




It follows that 

^. L in ss(TT^J ■ 

|X+|=n 

Below we shall use this formula from the right to the left in order to rewrite sums over 
partitions as multiple integrals. 

In order to apply (105) we split the sum in (99) as 

m 

E =£ E E • (106) 

({X+},{X-})ep2{M «=0 {mA^-})^P2m ({^+},{^-})ep2W 
({^+},{r})eP2© n=n \k+\=n 

%+\+\K-\=Nl2 

We would like to apply (105) to the rightmost sum in (106) after inserting it into 
(99). However, y„({A.+}|{^+})Z„({X+}|{^+}|{^-}) is not analytic in the variables 
X+ inside C", but for every has simple poles at = for k = 1, . . . ,n (the 
function G(?i,^) has a simple pole at ?i = ^ (see the appendix), the simple poles of 
detM(Xj',^^) at = are compensated by the simple zeros of b{X^)). Recall that 
the are free parameters. We choose them distinct from the Bethe roots Xj and all 
inside C. Then there is a simple closed contour F encircling all the but none of the 
Bethe roots Xj, and y„({X+}|{^+})Z„({X+}|{^+}|{^-}) is analytic inside C -T (see 
figure 2). Thus, (105) applies to C — F, and we obtain 
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E 

{{X+},{^-})ep2{X} 
\X+\=n 



Ynim\m)Zn{{i+}\m\m) 



n 



n 



{c-rr IfJi 2ra(l + a((0;)) J (l + a{^j )) J 

i'«({«;}}=i|{r})Z„({(0,}}=i|{^+}|{r}) • (107) 



This is already a multiple integral representation. Note, however, that it is still inap- 
propriate for performing the homogeneous limit and the Trotter limit for two reasons. 
First, the integrand still contains the auxiliary function ci(^J) which cannot be removed 
by applying a formula like (105) to the variables But the Trotter limit and the limit 
^ ^ in a(^) do not commute. Second, since the origin is a limit point of the sequence 
of Bethe numbers that determine the leading eigenvalue of the quantum transfer matrix 
in the Trotter limit, we could not avoid that Bethe roots would cross the contour T if 
it would include the origin which in turn would be necessary in order to perform the 
homogeneous limit. The way out of these problems is amazingly straightforward. In 
the next two sections we calculate the F-integrals explicitly and sum up the resulting 
terms. In this procedure all unpleasant terms ci(^j) cancel each other. 



3.3 Evaluation of the F-integrals 

Due to the symmetry in the COy of the integrand in (107) we may replace the integral 
with 



/. n n 



n-k 



n—k „ k 



7=1 



7=1 



(108) 



Our next task is to evaluate the F-integrals. This can be achieved by means of the 
residue theorem. We obtain the following result. 



f r-p-r d&n-k+j 

7r*l-}j27ii(l-Fo(co„_; 



-k+j)) 



= k\ Y, 

({r+},{4+-})€p2{^+} 
\^+-\=k 



y„({a),}}=il{^+})z„({co,-};;=il{^+}|{r}) 

. k 1 



)) 



n(^l+e a(^jn^^^,,^^_^^^^_^^^^ 



■ (109) 
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3.4 Resummation of the C -integrals 



Inserting (107)-(109) into (99) we obtain an expression that contains only integrals 
over the canonical contours C. Yet the summation with respect to the inhomogeneity 
parameters looks more complicated than before. It turns out that these sums can be 
simplified considerably by applying the following 

Lemma 5. For any function F(^) the identity 



k=0 (K+},{^-})ep2{^} 



n(l+eM^7)F(^7))l [n(l-e^F(^+)) 

7=1 y=i 



1^1 



7=1 



holds. 

Thanks to this lemma we end up with the equation 

1 



m „{m—n)((> 

M9|{^})=E^ I 



jJi 2jti(l + a(a)^)) 



n 



sh(a)^--^+-ri) 



detM(coy,^+) detG(a)y,^+) 

(111) 



for the finite Trotter number approximant to the generating function (35) of the a^-a^ 
correlation functions. Here we introduced the new abbreviations 



(112) 



M(o),,^, ) = t{^, '«i)n,h(eo,-to,-Ti)+^^^^'^^ nsh(to,-to, + Ti)- ^''^^ 

The integral on the right hand side of (1 1 1) is analytic and symmetric in the in- 
side C". It vanishes if any two of the are identical. The function l/b(A,) is meromor- 
phic inside C with only simple poles located at X = ^j. Hence, we can apply (105) to 
any simple closed contour F which lies inside C and encircles all the inhomogeneities 
^y. It follows that 



m (m— «)(p 



sh(cOy-^fe-Tl) 



n 



27ii(l + o(a)y))_ 



This is our final formula for the finite Trotter number approximant 4>A^((p| {^}) . We are 
now in a position to take the homogeneous limit and the Trotter limit analytically. 
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3.5 Taking the limits 



The information about the inhomogeneities is contained in the function b(k). From its 
definition we see immediately that 



lim b(X) 



sh(;i) 

sh(X-ri) 



(115) 



The information about the discreteness of the Trotter decomposition, on the other hand, 
is implicit in the function o(A,). Taking the Trotter limit means to take a(X) as a solution 
of equation (63). Thus, we conclude and summarize our main result. 

Theorem 1. The generating function of the cr^-o^ correlation functions of the XXZ 
chain has the following integral representation: 



(exp{(pL;r=l^ni}>r/. 



m g(m— n)<p 



n=0 



.|.Vr2m 



d^j [ s h(C;-Ti) 
sh(C;) 
^ sh(a)y-C;fc-Tl) 



sh(tOy) 



c 27ii(l + a(toy)) Vsh(coy -T|) 

detM(a)^-, ^k) detG(co^-, Cfe) , (1 16) 



where a(X) is the auxiliary function (63) that determines the free energy in the Trotter 
limit, 

and the function G{X, Q is the solution of the linear integral equation 

d(0 sh(2ri) Gioi.Q 



(117) 



G(U) = .(U)+/^ 



c 2%i sh(A,-to + r|)sh(A,-to-T|) l + a(to) 



(118) 



The contour C is the canonical contour in the non-linear integral equation formalism 
{see figure 1). Y is any simple closed contour which lies inside C and encircles the 
origin. 

As we have seen in sections 2.7, 2.8 both auxiliary function, a and a, can be equiv- 
alently used to express physical quantities. To switch from one representation to the 
other one can use the identity 



1 



l + a(;i) 
Inserting it into (118) and using 

sh(2Ti) 



1 



1 



l + a(;i) 



sh(X-C + Tl)sh(X-C-Tl) 



= r(X,Q+r(C,X) 



(119) 



(120) 
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we find, for instance. 



sh(2Ti) 



G(to,Q 



c 2%i sh(A,-to + T|)sh(A,-to-T|) l + a(to) 



(121) 



which equally well defines the inhomogeneous density function G(X, Q. The task of 
re-expressing the multiple integral representation (116) in terms of a is, of course, 
more cumbersome. One may, for instance, start from (114), insert (119), evaluate the 
integrals which do not contain a, and resum the resulting terms in a similar way as in 
section 3.4. One arrives at the following 

Corollary 1. When expressed in terms of a the generating function takes the form 



(-1)" 



«=0 



(n!) 



1^2 



n 

■j=i 



J^,^sh(g+ri) 
27ii 



Jco, 



n 



sh(a);-C^ + ri) 



sh((0^-) 



2;i:i(l + a(a)^-)) Vsh((Oy+ri)_ 
detM(tO^-, detG(tO^-, , (122) 



where M{Gij, C,k) and G{G)j, C,k) are the same as in theorem 1. 

In the remaining sections we shall discuss several special cases and limits of (1 16) 
and (122). 



3.6 One-point function 

We begin with the one-point function obtained from (116) or (122), respectively, for 

m = 1. In this case the sums on the right hand side of (1 16) and (122) consist of only 
two terms, the first of which does not contain any integral. In the second terms the 
integration over Y can easily be performed, and one remains with the formulae 

J(0 G((o) _ , , „(p, /■ d()i G((o) ^^^^^ 



UoeWi^r,;, e ^ Vc27ril + a(a)) ^ ^ Jc 



2ni 1 + a((o) 27ii 1 + o(a)) 

Here we inserted the identity G(co) = G(co, 0) which is obtained by comparing (72) 
and (118). 

Setting (p = in the above equation both expressions for the one-point function 
yield 1 as it has to be. Setting (p = in, on the other hand, and multiplying (123) by j 
we obtain two formulae for the magnetization, 

/^A _ 1 f d(o G(co) _ 1 f d(o G(a)) 

\^^T,h 2 ic27iil + o(to) 2^yc2ml + a(to)- ^ ^ 

These formulae are the same as obtained in section 2.8 based on a completely differ- 
ent reasoning, namely by taking the derivative of the free energy with respect to the 
magnetic field. 
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3.7 Emptiness formation probability 

In the limiting cases cp — ±00 (for cp — 00 one has to multiply by e^'^f first in order to 
obtain a finite result) our formulae (116), (122) simplify considerably. The resulting 
expressions on the left hand side of (116), (122), (ei | . . . em\)T,h or (^12 • • • ^m^T,h^ re- 
spectively, describe the probability to find a string of up- or down-spins of length m. 
For spin chains these probabilities where introduced in [13]. They come under the 
name of the 'emptiness formation probability' and recently became popular, since the 
integrals in the zero temperature multiple integral representation of the emptiness for- 
mation probability [4, 13] could be evaluated explicitly for small m [1]. 

For finite temperatures our formulae (116) and (122) yield the two alternative ex- 
pressions 



{e\\...em\) 



T,h 



(m!) 



1^2 



JJiJr 2ni V sh(g) J Jc 2ni{l + o(coy)) Vsh(cO;- -r]) 



sh((o^- 



yPi sh(g - Cfc - n) sh(0); -03k- T\) 



m 

= 1 

n=0 



(-1)" 

(n!)2 



.|=Vr27ii 



det ? (Hj) detG(tO^-, (125a) 
d(Sij ( sh(a)y) 



./Jl - Cfe + n) Sh(t0y - tOfe - Ti) 



c 27ii(l + o(to^-)) Vsh(to^ +T|) 
detf(Cfe,co^) detG(co^-,Cfe) (125b) 



for (ei I . . . Cm } ) y ^ and, similarly. 



(_l)m 



(m!) 



1^2 



^CjY^h(Cj+ri) 



n 



n^ir27tiV sh(Cy) y ic27ti(l + o(tOy)) Vsh(tOy + Ti) 
sh2(tO;-C;fe + Tl) 



sh(a)j' 



= 1 



A^ij sh(Cy - Cfe + n) sh(tOy - tOfe + Ti) 
1 

)2 

I-; 



det?(0)^-,^^) detG((0;,C^) (126a) 



f jg/sh(g-Tl) 



^ Jr 2ra \ sh(^ 



sh((jOy) 



■ " sh((D;-^fc + ri)sh((D — ^;^-ri) 
./fci 1 sh(C; - + 11) sh(a)^- - (0;t - Ti) 



c 27ii(l + a((0;)) \sh((Oy-ri) 
det ? (tO;, X^k) detG(tOy, • (126b) 
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3.8 Zero temperature limit 

How can we make contact with the zero temperature results obtained in [5]? It is 
clear by inspection that the auxiliary functions a and a defined by (63) and (65) do not 
approach a limit as T goes to zero. The inhomogeneities are unbounded in T. In order 
to obtain a sensible result we have to multiply by T first. Setting 8(A,) = — rin(a(A,)) = 
rin(a(?i)) the non-linear integral equations turn into 



2/sh2(Ti) r da sh(2Ti)nn(l + e-^('")/^) 



2ysh"(r|) r da 



sh(A,)sh(A,+T|) Jc 27ti sh(A,-Q) + r|)sh(A,-a)-ri) 



(127a) 



27sh2(n) f d(D sh(2ri)rin(l+e^N/^) 
sh(A,) sh(A, - Tj) Jc 2ni sh(A, - co + rj) sh(A, - to - r|) ' 

Here we shift the left and right edges of the integration contour to ±T|/2, where 
e{X) is real (recall that the contributions from the upper and the lower edge cancel 
each other). We see that it only matters if e(k) is positive or negative as T goes to 
zero. In (127a) only the negative part of e{X) on the contour contributes to the zero 
temperature limit, in (127b) only the positive part. This divides the contour C into two 
disjoint pieces C^'^^ and C^~^ on which e(A,) can be calculated by means of the linear 
integral equations 

Em = h+^^^^^- f sh(2Ti)e(a)) 

sh(;i)sh(;i + ri) Jci-)2nish{l-(o + T])sh{X-(o-'(])' ^ ^ 

cm-!, 2/sh^(Ti) f da sh(2Ti)e(co) 

sh(X)sh(A.-T|) 27ii sh(X-to + Ti)sh(X-to-Ti) ■ 

These integral equations can be transformed into integral equations for dressed ener- 
gies, which is, however, not our actual aim. The most important conclusion instead is 
that 

1 1 r-^o fo onc(+) 



l + o(;i) l+e<^)/T ' ]l on C^-y 



(129) 



This function behaves like a 'Fermi function on c:'. 

Applying (129) to the integral equation (121) and shifting the arguments appropri- 
ately we obtain 

^ 2'^ 2J ^(-)+E27tish(X-to + 'n)sh(A.-to-'n) ^ ^ ^ 

This equation has to be compared with equation (2. 16) of [5]. The two equations agree 
if we set 

G(X-n ^-n)=27tip(X,Q (131) 

and identify the integration contour in equation (2.16) of [5] with — (c*^^-* + ^). A 
close inspection (for details see e.g. [10]) of (127), (128) shows indeed that (for h>0. 
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T| < 0!) the contour is a line segment ReA, = x\/l, —a < ImA, < a (where a > 
depends on the magnetic field h). This is sufficient to identify (130) with (2.16) of [5]. 

Let us introduce the shorthand notation I = —{C^^'' + Then, using (129) and 
(131), our second expression (122) for the generating function of the o^-o^ correlation 
functions can be evaluated in the zero temperature limit. We obtain 

lmi(exp{(pL-^ie„i}>y^ 



1 



n2 



j^Vr+5 27riVsh(C,-n); A ^Vsh((0,- + ^), 

detM(tO^-,C;fe) detp(tOy,C;fe), (132) 



^ sh(a)^--^^ + ri) 



where p(A,, Q is the solution of the linear integral equation 

sh(2ri)p((o,g 



-27ti 



dm 



sh(?i - CO + ri) sh(X - CO - ri) 



(133) 



which is precisely the result of [5]. 



3.9 High temperature limit 

In the limit of infinite temperature all correlations should vanish. How can this be seen 
from our formulae? For T — > oo the inhomogeneity on the right hand side of the integral 
equations (63) and (65) vanishes. The remaining homogeneous equations are solved 
by a{X) = 1 or a{X) — 1, respectively. Inserting this, for instance, into (116) we see 
that the only singularities in the resulting integrand are the simple poles of the function 
G(co, Q. Therefore the integrals can be calculated by means of the residue theorem. 
Using similar techniques as described in section 3.4 we can simplify the result which 
finally turns into 

hjn(exp{cpi:-=ie„i}>^,= (^^j . (134) 
Applying (36) we obtain 

lim(afc4>j,, = (135) 

as we expected. 



4 Conclusions 

In this work integral representations for finite temperature correlation functions of an 
integrable model have been derived for the first time. We concentrated on the gen- 
erating function (35) of the o^-a^ correlation functions, since it appeared particularly 



30 





Im 
2 










r, 




c 






j 

c 












2 



Re 



-y 



Figure 3: Contours of integration in the critical regime for < A < 1. For —1 < A < 
one lias to choose |y| < | — Alternatively one may switch the sign of / making use 
of the equivalence of the parameter regimes —1 <A<0, y>0 and 0<A< 1,7<0. 



appropriate to us for developing the method. Now this method may be applied to other 
generating functions and also directly to two-point functions. Work in this direction is 
on the way. 

It will be very interesting to further study the integral representations. How can 
they be utilized to obtain the long-distance asymptotics? How efficient are they for 

numerical calculations? Are there special cases where part of the integrations can be 
performed analytically? We hope to address these questions in forthcoming publica- 
tions. 

We restricted our presentation to the off-critical regime A > 1. We would like to 
stress, however, that our results straightforwardly extend to the critical regime |A| < 1, 
where only the canonical contour has to be redefined (see figure 3). The situation is 
much the same as in the zero temperature limit [5]. 

Acknowledgement. The authors would like to thank M. Karbach for helpful dis- 
cussions. 
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Appendix: A ratio of determinants 

Here we prove lemma 3. We define a column vector n(^) with entries 

nji?,) = t{ljX)-t{?,,lMt,), j=l,...,N/2 (A.1) 
and another column vector 

m=No'ni^). (A.2) 

The matrices A'^o and A'^„ have the last A''/2 — n columns in common. Thus, we conclude 
that 

= dQtNQ^Nn = det(j(|i), . . . ,j(|„),e„+i, . . . ,e;v/2) = det (ey, , (A.3) 

where the ej are the canonical unit column vectors having a single non-zero entry 1 
in the jth row, and j and k on the right hand side run from 1 to n. In (A.3) we have 
expressed the ratio of two N/2 x N/2 determinants as a single nxn determinant, a 
trick we borrowed from [6]. 

In the second step of our proof we transform the equation A^oj(^) = n(^) that de- 
termines j into an integral equation. For this purpose we first of all rewrite it in com- 
ponents, 

-a'{h)h{^) + L^l^^r^T^K^, ; = t{hX)-t{^MH^) . (A.4) 

iti sh{lk - ^/ + Tl) sh{lk -h-n) 

We define a function 

/=i sh(A, -Xi + r\) sh(A, - A,/ - r|) 
and observe that this function determines j (^) through 

F{Xk,^,) = a'{h)m)- (A.6) 

We shall now assume that ^ is located inside our canonical contour C sketched in 
figure 1. Then the only singularity of F(X,^) as a function of X inside the rectangle Q 
(see again figure 1) is a simple pole atX = b,. The residue at this pole is 

res{F(X,^)}|^ . = limsh(X-^)F(X,^) = -1 -a(^). (A.7) 

This residue gives an additional contribution when we apply equation (105) to (A. 5) 
in order to transform the sum on the right hand side into an integral. We end up with 



32 



The latter equation suggest to introduce a function 



l + a(^) 



(A.9) 



This function then solves the linear integral equation (97) defined in lemma 3. Coming 
back to our original notation (see (93)) we may finally express the ratio of the two 
determinants as 



detNn 
detA/o 



n 

■1=1 



i+a(^r) 

a'(V) 



detG(^+^+), 



(A. 10) 



and our lemma is proven. 

Note that is follows from (A.7) and (A.9) that the only singularity of G(k, ^) (con- 
sidered as a function of X) inside the rectangle 2 is a simple pole with residue — 1 at 
"k — t,. This fact is frequently used in the derivation of the multiple integral represen- 
tations (116) and (122). 
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